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1. Introduction 

Precision measurements of the CKM matrix put the Standard Model to a stringent test and 
constrain possible physics beyond it. Using the measured frequency of B q — B q , q € {d,s} oscil- 
lations to determine the CKM matrix elements \V tq \ requires a reliable lattice calculation of the 
non-perturbative B q — B q mixing matrix elements \m\ fg Bg q - A 2 + 1 flavor, unquenched calcu- 
lation of fs and Bb has been carried out by the RBC-UKQCD collaboration in the infinite heavy 
quark mass limit using light domain- wall fermions on a (2 fm) 3 spatial volume [1, 2]; this is cur- 
rently being extended to a (3 fm) 3 spatial volume and towards physical light quark masses [3], In 
the following, we discuss the perturbative lattice-continuum matching of the operators relevant for 
the RBC-UKQCD calculation, following in part the detailed discussion in Refs. [2, 4]. We also 
point out the subtle degeneracy of heavy-light meson ground states, and discuss its implications for 
the extraction of /# and Bb from lattice correlation functions. 



2. Action and Feynman Rules 

The heavy b quark is described by an improved lattice version of the static limit of heavy quark 
effective theory with smeared, SU(3)-projected gauge links Vo(x,t) to reduce noise: 



^ h(x,t + a) h(x,t + a) — Vo(x,t)h(x,t) 



(2.1) 



The SU(3) projection (discussed in Ref. [2]) simplifies perturbative calculations by allowing the 
smeared gauge links to be expanded in terms of an effective gauge field B c '(x,t); in momentum 
space By(q) = h^(q)A"(q), where A°(q) is the physical gauge field and h^(q) is a form factor 
depending on the smearing scheme. We focus on one of the two schemes used in the RBC-UKQCD 
calculation (one-level APE blocking with parameter a = 1), resulting in a heavy quark gluon vertex 

= -igoT a d^e- iiko+ ^ /2 Fj(*,tf) = -igofh^e-to+W 2 , (2.2) 
where go is the bare lattice coupling, q is the gluon momentum, and h^{q) is given by 



= (h (q), hj(q)) = ^1 - \ t Sin 2 (f ) ,\ sin (f ) sin (|j j . (2.3) 

The heavy quark two-gluon vertex and the heavy quark propagator are given in Ref. [4]. 

The light quarks are described by the domain-wall fermion action. Each light flavor is repre- 
sented by a (4+ 1) -dimensional Wilson-style fermion field \j/ s (x,t) where 1 < s < N labels the 
coordinate in the fifth dimension. The physical quark field q(x,t) is constructed from chiral surface 
states at s = 1 and s = N via q(x,t) = PR\y\(x,t) +PLYN(x,t). The domain-wall height M$ is a 
fixed parameter of the theory; we set M=, = 1.8 to match the RBC-UKQCD calculation. A detailed 
description of domain-wall fermions and their perturbative treatment for our choice of gauge action 
is given in Ref. [4] and references therein, especially Ref. [5]. In the perturbative calculation the 
light quark masses were set to zero and the size ,/V of the fifth dimension was taken to be large, 
resulting in an exact chiral symmetry as ,/V — > °°. The gluons were described by the Iwasaki gauge 
action, whose Feynman rules are given in Ref. [4]. 
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3. Perturbative Lattice- Continuum Matching at One-Loop 

The full QCD operators relevant for the extraction of fs and Bb, defined in MS(NDR) at the 
scale Hb = mb of the b quark mass, are the axial vector current A p = by p y^q and the parity-even 
part of the AB = 2 vector-axial four-quark operator: 

[byP{\ - y 5 )q] [by p {\ - y 5 )q] -> O vv+ aa = {by p q) {by p q) + {by p y 5 q) {by p y 5 q) . (3.1) 

We match these operators at the scale \lb to lattice operators in the static effective theory (de- 
scribed in Sec. 2) at the lattice scale aT x via the continuum version of the static effective theory 
renormalized at a scale fx. Throughout our one-loop calculation we choose to set pi = a~ l ; in the 
RBC-UKQCD calculation, the lattice scale is given by a- 1 = 1.62 GeV. The full QCD operators 
are related to continuum static operators by 

A p (to) = C A (to,p)Ap(ji) + 0(/iqcD/lib), (3-2) 
O vv+AA {pLb) = Zi(Hb,LL)d vv+AA {iJ.) +Z2(Hb,n)d ss+P p(n) + ^(AgcD/Atfo)- (3.3) 
In terms of the static quark and antiquark fields h^{x) = e±»"* v *(i ± f)b{x) jl and for ni/, — > oo, 

A p =h i+) y p y 5 q, (3.4) 
5 VV+AA = 2 (h [+) y p q) (h H y p qJ +l{h [+) yPy 5 q) (h { - ] y p y 5 q) , (3.5) 

Oss+PP = 2 (h i+) q) (h^q) +2{h [+) y 5 q) (h H y 5 q) . (3.6) 

The static effective action discussed in Sec. 2 describes M + ) with v = (1,0), corresponding to 
a stationary meson. The constants Ca(jU/,,ju) and Zi^jU^jU) are known at one-loop; they are 
summarized in Ref. [2] . Using the latest PDG values for a^ s (mz) and nib, and running the coupling 
down at four-loops with the physical number of flavors to determine a 4 MS (ji^) and a^ s (/x) we 
obtain C A = 1.057, Z x =0.934, Z 2 = -0.151. 

We now describe the matching A p (\i) = Ca(h ,a~ l )a~ 3 A l p t of the heavy-light axial currents 
Ap(ju) and A 1 ^ 1 (which is dimensionless) in the continuum and lattice versions of the static ef- 
fective theory. Results for the four-fermion operators are summarized at the end of this section. We 
compare the correlation function ({h{x)Tq{x))h(y)q{z)) in both theories; in this discussion only 
one heavy quark field = h enters. For the axial current T = y p y$, but the light quark chiral 
symmetry and the heavy quark spin symmetry h — ► e^"^' E ' k,ak 'h of both the continuum and the lat- 
tice theory render the matching T-independent. At one-loop and for small external quark momenta 
p ~ the continuum and lattice correlation functions are 

{{ h Yq)h q) = ^r(i + 8v)% ((hr q )hq) M = g^r(i + gO (1 " w °f H ' Z ^ , (3.7) 

ipo if ipo if 

where the Feynman diagrams contributing at one-loop are shown in Fig. 1. All Z- factors have 
values 1 + &{a s ), and the vertex corrections SV, SV lat are &{a s ) and T-independent as noted above. 
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(c) 



The continuum quantities are known [4]; we focus 
on the lattice correlation function: wo = 1 — M$ is a 
domain-wall fermion specific constant, and an over- 
lap factor 1 — Wq connecting the five-dimensional and 
physical quark fields is present even at tree level. The 
light quark wavefunction renormalization Z^Zj due to 
Fig. 1 (a) and (b) was calculated in Ref. [5]. zJf 1 can 
be viewed as the four-dimensional wavefunction renor- 
malization, while Z w renormalizes the overlap factor 
1 — Wq. Due to tadpoles, the one-loop correction to Z w 
is enormous. As described in Ref. [5], this is remedied 
by reorganizing the perturbation series according to the 
mean-field approach, resulting in the prescriptions M5 — > M5 = M5 — 4( 1 — u) , wq 




(d) 




(b) 



(e) 



Figure 1: One-Loop Corrections to the 
Heavy-Light Axial Current 
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i/2^iat to k e ma( j e throughout the calculation; here u = P 1 / 4 where P is 




the measured average plaquette (for the RBC-UKQCD calculation u = 0.8757) and the superscript 
'MF' identifies mean-field improved quantities. We calculate the matching factor Ca^cT 1 ) using 

both the usual continuum MS coupling and a mean-field 
improved version, enabling an estimate of &(ctj) correc- 
tions. a^ s (/i) was obtained by running down to the c 
quark mass with the physical number of flavors and back 
up to ji, using only three dynamical flavors to match the 
RBC-UKQCD 2 + 1 flavor calculation: af S {pi) = 0.326 
and a^ F (^) = 0.177. The calculation of the vertex cor- 
rection 8V M in Fig. 1 (c) and the heavy quark wavefunc- 
tion renormalization Zf l in Fig. 1 (d) and (e) is straight- 
forward [2, 4]. Infrared divergences only occur in QED- 
like diagrams and are regulated by a gluon mass X which 
cancels from the matching factor. Furthermore, only the 
unsmeared part of h^(q) in Eq. (2.3) gives rise to infrared divergences; the sine functions in 
the smeared part of h^{q) cancel all infrared divergent loop propagators. A generic feature of 
domain-wall fermion perturbation theory is the appearance of correlation functions (q(— p)y s (p)) , 
(Ws{—p)q~{p)} connecting external four-dimensional quarks to five-dimensional quarks propagat- 
ing in loops, as shown in Fig. 2. A subtlety pointed out in Ref. [6] is that the correct renormalization 
prescription for z^ atMF includes the linearly divergent heavy quark mass renormalization: 



Figure 2: One-Loop Vertex Correc 
tion to the Heavy-Light Axial Current 



7 lat, MF 



1 



.dL(po) 



dpo 



po=0 



(3.8) 



where the heavy quark self energy £(/?o) itself is not affected by mean-field improvement. Com- 
paring the correlation functions in Eq. (3.7) after mean-field improvement gives a matching factor 



C A (pL,a l ) 



MFj2)zMF 



1 + ^(-1.584). (3.9) 
3K 
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The overall factor Z<s,{jjLb,a~ l ) = Q^^Ca^gT 1 ) relating the axial currents in full QCD and 
the lattice static effective theory, computed using both a 4 MS (jU) and a^ ¥ {pl), is Z^ s (/Xfo,a _1 ) = 
0.902, Z^ F (/i/,,c? _1 ) = 0.961. While the one-loop result is small and reliable, the large difference 
between a] VIS (/i) and a^ ¥ {pl) induces a ~ 7% systematic error ultimately warranting nonpertur- 
bative renormalization. 

For completeness, we quote the lattice-continuum matching constants (calculated in Refs. [2, 4]) 
for the four-fermion operators in Eqs. (3.5) and (3.6). For i £ {VV +AA,SS + PP} and at one-loop 

= (1-( w mf ) 2 )z mf Z - MF ^^' 1 )^ 6 0!" Z^ +AA = l + ||(-4.462), Z™ PP =l, 

(3.10) 

where the C>J at are dimensionless. Since the coefficient Z2 of Oss+pp in Eq. (3.3) is &(a s ), only 
the domain-wall overlap factors contribute to the lattice-continuum matching for this operator. 
While formally inconsistent, we use the one-loop mean-field improved values of the overlap factors 
throughout to ensure tadpole-safety. Combining Eqs. (3.3) and (3.10) we get: 

Ovv+aa = ZvA(lJ-b,a l )a 6 Oy V+AA +Zsp(plb,a l )a 6 O l ^ +PP , (3.11) 
Z^f= 0.902, Zfl = 0.769, Z^?= -0.123, Zf P F = -0.133. (3.12) 

4. Ground State Degeneracies of Static-Light Mesons and fg, Bb on the Lattice 

Let H be the Hamiltonian corresponding to the full lattice action in Sec. 2. For any t, the heavy 
quark action in Eq. (2.1) is invariant under h(x) — > e' e ^h(x) for a set of V /a 3 parameters 6(x), 
where V = L 3 is the spatial lattice volume. If &(x) is the generator corresponding to 6{x) then 

[e(x),h(y)] = h(x)8^ [e(x),h(y)]=-Ht)Sxy, [©(*)> ®60] =o, [©(*),#] = o. (4.i) 

Simultaneously diagonalize H and all &(x). Since Eq. (4.1) implies that h(x) and h(x) raise 
and lower the eigenvalues of &(x) by 1, and the charge conjugation invariance of QCD implies 
&(x) |0) = 0, the spectrum of &(x) contains Z. Define the unit-norm state to be the lowest 

energy state with the quantum numbers of a B meson which also satisfies &(y)\B(x)) = 5^t\B(x)}. 
Thus (B(x)\B(y)} = 5^, and we can interpret these states as having the heavy quark localized at 
a fixed lattice site with the light quark smeared out around it. Since T(i)&(x)T(i)~ l = ®(x+i), 
where T(i) is a lattice translation by a in the spatial direction i, all B meson ground states \B(x)} are 
degenerate. We also define total spatial momentum eigenstates |B(&/)), where G Z (z = 1,2,3): 

\B(h)) = V^£e- rk '-*\B(x)), Z, = ^(Zi,fe,Z 3 ), -Ya <k -ka' ^(k v )\B(k)) = 2V8 n . 

(4.2) 

As a — > 0,V — > 00, these states reduce to continuum momentum eigenstates \B(p)) c with conven- 
tional static effective theory normalization c (B(p')\B(p)) c = 2{2k) 3 8^\p' — p). In the mi, — > 00 
limit, these states only differ from the corresponding full QCD states by a factor of ^/rag. Thus: 
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fsV^B = <0|A (0,0)|fl(p = 0)) c = Z^ F a- 3 (0^(0,0) \^"L \ B W) 

= V2Z™ F a- 3 / 2 (0\A%%0)\B(d)} = ^Z^aT 3 ' 2 ^. 

(4.3) 

In complete analogy to the above, we can construct B meson ground states states \B(x)). Using 
these and Eq. (3. 11), the calculation of the B — B mixing matrix element I m B f B Bb = (B\0w+aa\B) 
is reduced to the calculation of the lattice quantities (S(0)|O| at (0,0)|B(0)), i € {VV +AA,SS+PP}. 

The degeneracy of the states \B(x)) complicates the extraction of Og 1 and (fi(0)|O- at (0)|fi(0)), 
since even a large time separation of source and sink may not project onto a unique B meson ground 
state: different combinations of the \B(x)} may enter the correlation functions used for calculating 
the matrix elements and those used for normalization. To see this, consider the extraction of <J>£ £ ; 
we now work exclusively in the lattice theory. Define local and smeared B meson interpolation 
operators A%(x,t) = h(x,t)y Q y 5 q(x,t), A s Q (t) = ZfeAvLteAvHyrfYoYslfct), where AV is a fixed 
subvolume of V and the smeared operators are Coulomb gauge fixed. From experience, local- 
local correlation functions in the static effective theory are prohibitively noisy; instead calculate 
the local-smeared and smeared-smeared correlation functions. Inserting a complete set of states 
\B(w))(B(w) | + (higher energy states) with the correct quantum numbers, we have as t — > °°: 




^(f) ^^(0\A L (lt)A s (0^\0)=^fe- m s' £ (B(w)\ £ £ ?(?,0))t>)5*(Z,0)|0> , (4.4) 



xev 



veAVzeAV 



V ss (t) = (0\A*(t)A> (0y\0) = e- 



,weV 



I (B(w)\ £ £ q(y,0)TbYsh(z,0)\0) 

yeAVzeAV 



(4.5) 



where m* B is the unphysical mass of the lattice B meson. Since c to SS (t) contains a sum over squares, 
the use of a naive ratio ~ ^ LS {t) /^ ss (t) requires a translationally invariant wall source AV = V 
to project onto the unique state of zero-momentum. In this case the sums over w only give a 
factor of V/a 3 and <J>jf = tf LS (t) / yV^f )«r m * f V/a 3 . To remedy the poor overlap of the wall 
source with the B meson ground state - especially on large lattices - consider a fixed box source 
and a series of box sinks summed over an entire timeslice to project onto zero momentum; this 
approach also allows more general types of smearing, such as the use of an atomic wavefunction. 
Let A s Q {w,t) = EyeAV-YtzeAV-Hy^YoYsqiZit) where AV$ is a box of fixed size located at w and 
AVg = AV, Aq(0,?) =Aq(?). Define a corresponding smeared-smeared correlation function and 
insert a complete set of momentum eigenstates 977 £r \B(ki))(B(ki) \ + (higher energy states); then 



as t 



^(0\A s (w,t)A s (0y\0) 



2V 



-£(0|Ag(w,/)K0))(B(0)^(0)t|0) 



2a 3 



<a(o)|Ag(o)t|o) 

(4.6) 
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Since \B(0)) = (2a 3 ) 1 / 2 ^ vi ; \B(w)), we can rewrite the right side of Eq. (4.4) and obtain another 
ratio for <J>g l which reaches a plateau more quickly due to the improved ground state overlap: 



cgis^mitli I ^ SS{t) = ^ { B(d)\A s (0Y\0) /y |<S(0)|Ag(0)t|0>| 2 = (4.7) 
The calculation of (B(0)|O| at (0)|S(0)), i G {VV +AA,SS + PP} is considerably simpler. Define 

V ,(T,t) = £(0\A S (T)Ol*\x,t)A s (0y\0), (4.8) 

xeV 

where A (T) =Y,veAvllzeAvq{yiT)YoY5h(z,T). Proceeding as above, we have as t, T — t— > oo\ 



(B(0)|Ol at (0,0)|B(0)) =tf 0l (T,t) /^ SS {T) = c g 0i {T,t)e m * TI2 ^tf ss {T -t)tf ss {t). (4.9) 

Here no zero momentum projection is necessary; the use of c to SS for smaller time separations simply 
reduces noise. Using Eqs. (4.7) and (4.9) we can thus calculate fg and Bb using only box sources 
and sinks. These are preferable to wall sources, whose poor ground state overlap led to late plateaus 
in the V = (2 fm) 3 RBC-UKQCD calculation and presents an even bigger problem for the ongoing 
extension to V = (3 fm) 3 . It is worth emphasizing that this simple method relies on the particular 
properties of the static effective theory, and further such improvements might be possible. 
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